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Abstract 

We study a Gaussian relay network in which multiple source-destination (S-D) pairs communicate through 
relays without direct links between the sources and the destinations. We observe that the time- varying nature of 
wireless channels or fading can be used to mitigate the interference. The proposed block Markov amplify-and- 
forward relaying scheme exploits such channel variations and works for a wide class of channel distributions 
including Rayleigh fading. We completely characterize the degrees of freedom (DoF) region of the Gaussian relav 
network. Specifically, the DoF region of the if-user M-hop Gaussian relay network with K ra nodes in the m 
layer is the set of all (di, • • • , dx) such that d{ < 1 for all i and J2i=i d% < min{ifi, • • • , Km+i}, where di is 
the DoF of the z" 1 S-D pair and K — K\ = Km+i is the number of S-D pairs. We further characterize the DoF 
region of the Gaussian relay network with multi-antenna nodes and general message sets. The resulting DoF regions 
coincide with the DoF regions assuming perfect cooperation between the relays in each layer. 

O . I. Introduction 



Characterizing the capacity of Gaussian relay networks is one of the fundamental problems in network information 
theory. However, for Gaussian relay networks, the signal transmitted from a node will be heard by multiple nodes 
l/~> ■ (broadcast) and a node will receive the superposition of the signals transmitted from multiple nodes (interference) and 
there exist fading and noise, which make the problem complicated. To overcome such difficulties, simplified wireless 
network models have been developed in [1], [2], [3], [4], [5], [6] that provide intuition towards an approximate 
capacity characterization of single-source Gaussian relay networks [7]. 

Unlike the single-source case, the capacity or an approximate capacity characterization of multi-source Gaussian 

O ! relay networks is very challenging since the transmission of other sessions acts as the inter-user interference. Due 

to the interference, the extension from the result in [7] is not straightforward. Recently, remarkable progress has 

been made on multi-source problems in [8], [9], [10], [11], [12] and the references therein. It was proved in [9] 

that the Han-Kobayashi scheme indeed achieves the capacity of the two-user Gaussian interference channel within 

^h ' one bits/s/Hz. The capacity of the if -user Gaussian interference channel has been characterized in [10] as 

yl0g(P)+0(k>g(P)) (1) 

if channel coefficients are sufficiently independent and drawn from a continuous distribution, where P denotes the 
signal-to-noise ratio (SNR). To show the degrees of freedom (DoF) or capacity pre-log term of K/2, the technique 
of interference alignment was used, which minimizes the overall interference space by aligning multiple interfering 
signals from unintended sources at each destination. The concept of interference alignment has also been used 
to characterize the DoF of the if -user multi-antenna Gaussian interference channel [11] and the X-network in 
which each source has independent messages for all destinations [12]. Another alignment technique called ergodic 
interference alignment has been proposed in [13] showing that, for a broad class of channel distributions, half of 
the interference-free ergodic capacity is achievable for each user in the if -user Gaussian interference channel at any 
SNR. Exploiting the inseparability of parallel interference channels [14], [15], the ergodic interference alignment 
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Fig. 1. Example of a two-user two-hop Gaussian relay network. 



scheme jointly encodes messages over two specific channel instances to align the interference. A similar concept 
has been also applied for the finite field case in [13], [16]. 

Not only can the interference be aligned, but also it can be cancelled or partially cancelled for multi-hop 
Gaussian relay networks. Assuming amplify-and-forward (AF) relays, each destination may receive multiple copies 
of an interfering signal from different paths and potentially these copies can cancel each other through a suitable 
choice of the amplification factors of relays. Reference [17] has shown that partial interference cancellation using 
AF relays achieves the capacity of two-user two-hop Gaussian networks within a constant bit gap in some scenarios. 
Also, the interference can be completely removed so that the optimal DoF of K is achievable for if -user two-hop 
Gaussian networks if the number of relays is greater than or equal to K 2 [18]. 

In this paper, we study general multi-source multi-hop Gaussian relay networks. We observe that the time-varying 
nature of wireless channels or fading can be exploited to cancel the interference. As a simple example, consider a 
two-user two-hop Gaussian relay network in Fig. Q] in which 
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for even t, where H m [t] is the mr 1 hop channel matrix at time t. Note that if odd and even time slots are used 
separately, each S-D pair can only achieve 1/2 DoF since there is no path between the first S-D pair for even t 
and the second S-D pair for odd t. On the other hand, if the relays amplify and forward their signals with one 
symbol delay, then the interference can be completely cancelled since H2[t + l]Hi[t] becomes the identity matrix. 
Hence every S-D pair can achieve one DoF simultaneously. We generalize this idea to multi-source multi-hop 
Gaussian relay networks for a wide class of channel distributions including Rayleigh fading. The key ingredient is 
to set appropriate delays in AF relaying at each layer such that overall channel matrices become diagonal matrices 
with non-zero diagonal elements, which guarantees interference-free communication. Under this class of channel 
distributions, we completely characterize the DoF region of multi-source multi-hop Gaussian relay networks. This 
improves upon our previous result that showed total K DoF is achievable for if -user if -hop networks with K 
relays in each layer when K is even and a similar technique has been proposed for linear finite-field multi-hop 
networks (see the conference papers [18], [19]). We further characterize the DoF region of multi-source multi-hop 
Gaussian relay networks with multi-antenna nodes and general message sets. 

This paper is organized as follows. In Section [III we explain the underlying system model and define the DoF 
region. In Section [nil we state the main results of this paper, the DoF regions of Gaussian relay networks. In 
Section [iVl we propose a block Markov AF relaying scheme and derive its achievable rate region for two-hop and 
three-hop Gaussian relay networks. From the achievable rate regions, we characterize the DoF region of Gaussian 
relay networks in Section [V] We conclude this paper in Section [Vj and refer to the Appendix for the proof of the 
result in Section IIV-B1 

II. System Model 

In this section, we explain our network model and introduce encoding, relaying, and decoding functions. Based 
on this model, we define the capacity region and the DoF region. Throughout the paper, we will use A, a, and A 



TABLE I 

Summary of notations 



A T (ora T ) 


Transpose of A( or a) 


A^ora 1 ") 


Conjugate transpose of A( or a) 


||A|| F (or||a|| F ) 


Frobenius norm of A( or a) 


A" 1 


Inverse matrix of A 


[A]y 


(i,j) element of A 


[A], 


i row vector of A 


I 


Identity matrix 


diag(ai, • • • ,o n ) 


Diagonal matrix with [diag(ai, • • • , a n )]a = on 


real(a)( or imag(a)) 


Real (or imaginary) part of a 


|a| 


Absolute value of a 


card (A) 


Cardinality of A 



to denote a matrix, vector, and set, respectively. The notations used in the paper are summarized in Table H 



A. Gaussian Relay Networks 

We study a Gaussian relay network consisting of M + 1 layers with K m nodes in the rrv layer. We assume 
that the number of hops M is greater than or equal to two. The nodes in the first layer and the last layer are 
the sources and the destinations, respectively. Thus K = K\ = Km+i is the number of S-D pairs. We assume 
full-duplex relays so that all relays are able to transmit and receive simultaneously, but the results in this paper 
can be straightforwardly applicable for half-duplex relays by scheduling of hops. For notational simplicity, let us 
denote the v node in the rrv layer as node (z, m), where i £ {1, • • • , K m } and m £ {1, • • • , M + 1}. 

Consider the rrv^ hop transmission in which the nodes in the rrv^ layer transmit and the nodes in the (m + l) tJl 
layer receive. Let x^ m [t] denote the transmit signal of node (i, m) at time t and yj, m [t] denote the received signal 
of node (j, ra + 1) at time t. Then the input-output relation of the mr 1 hop is given by 



yj,m[t] 



£ 

2=1 



flji^m \P \ %i, m \p \ ~r Zj,m\P\i 



(4) 



where hj^ m \t] is the complex channel from node (i,ra) to node (j, m + 1) at time t and Zj jm [t] is the addictive 
noise of node (j, m + 1) at time t. We assume that Zj iTn [t\s are independent and identically distributed (i.i.d.) and 
drawn from A/c(0, 1) and each node satisfies the power constraint P. 

Let us denote x m [t] = [ai, m [*] 3 • • • , x K m ,m[t]] T and y m [t] = [yi, m [t], • • • ,yK rn + 1 ,m[t]] T , which are the K m x 1 
dimensional transmit signal vector and the K m +i x 1 dimensional received signal vector of the mr 1 hop, respectively. 
Then the rrv hop transmission can be represented as 



(5) 

i ,v ji,m\p\ 



y m [t] = H m [t]x m [t] + z m [t], 

where H m [t] is the K m+ i x K m dimensional complex channel matrix of the rmr 1 hop with [H m [t]]^ = h 
and z m [t] = [^i, m [t], • • • , ZK m+1 ,m[t]] T is the K m +i x 1 dimensional noise vector of the mr 1 hop. 

In this paper, we assume time-varying channels such that /iji, m [t]'s are i.i.d. drawn from a continuous distribution 
Ph(')- We consider a class of channel distributions such that Phip) > Phip) if M < H- The channel state information 
is assumed to be available at all nodes, i.e., each node knows Hi[£] to Hm[*] at time t. 



Ah 



B. Problem Statement 

Based on the network model, we define a set of length n block codes. Let W{ be the message of the z Ui source 
uniformly distributed over {1, • • • , 2 nR% }, where Ri is the rate of the i S-D pair. Then a (2 nRl , • • • , 2 uRk ; n) 
code consists of the following encoding, relaying, and decoding functions: 

• (Encoding) For i e {1, • • • , K}, the encoding function of the i^ source, or node (i, 1), is given by fi^t : 
{l,--- ,2 n ^}^C such that 

Xi,i[t] = /i,i, t (Wi) for tG {!,-•• ,n}. (6) 



• (Relaying) For m E {2, • • • , M} and i E {l,-" ^K m }, the relaying function of node (i,m) is given by 
fi,m,t ' C^ 1 -► C such that 

Xi,m[t] = /i,m,t(yi,m-l[l],'" >2/i,m-l[*- 1]) for t E {1, • • • ,7l}. (7) 

• (Decoding) For i E {1, • • • , if}, the decoding function of the z tJl destination, or node (i, M + 1), is given by 
# :C n "►{!.,••• ,2 ni *} such that 

Wi = 5i (Vi,M [1] , • ' • , 2/i,M N ) • (8) 

The probability of error at the z tJl destination is given by P e ^ = Pr(W^ 7^ W^). A rate tuple (i?i, • • • , i?#) is 
said to be achievable if there exists a sequence of (2 nRl , • • • , 2 ni?K ; n) codes with P e ^ — » as n —> 00 for all 
i E {1, • • • , if}. The capacity region C is the closure of the set of all achievable rate tuples. In the same manner 
as for the if -user interference channel [10], we define the DoF region as 



V= {(d u --- : d K )e 



»K 



VK--- ,wk)em^, 



K ( K R% \ \ 

y^Widi < limsup sup Y^Wi; — ^ f , (9) 

where d^ is the DoF of the r S-D pair. 

C. Multi-antenna and General Message Set 

We also study a more general case in which each node is equipped with multiple antennas and each source has 
messages to a subset of destinations. Let L^ m denote the number of antennas of node (i,ra) and W 9 denote a 
set of messages, where W g C {Wn, • • • , H / k m+1 k 1 }. Here, W^ is the message from the i^ source to the j^ 
destination and ifi ^ Km+i in general. Similar to Section HTBl the achievable rate region 1Z(W g ) and the capacity 
region C(W g ) can be defined. The DoF region is defined as 



V(W g ) = { {dji} Wji ew 9 e 



p card(W g ) 



ViwjijWjieWg € 



p card(W g ) 



V" Wjidji< limsup sup V" w J^r^B ] p ^ 10 ^ 

w ji ew 9 p ^°° \{^}^ievv 9 ec(w 9 ) WjieWg log^y j 

which is a simple extension of ©. Here, <i^ is the DoF from the r source to the j 1 destination. 

III. Main Results 

Throughout the paper, we characterize the DoF region of the Gaussian relay network. We simply state the main 
results here and derive them in the remainder of the paper. 

Theorem 1: The DoF region V of the Gaussian relay network is the set of all (di, • • • , dx) such that 

(U < 1 foriE {l,--- ,if}, (11) 

K 

J2d t < min {K m }. (12) 

f-^ mG{l,-,M+l} 

Proof: We refer to Section IV-AI for the proof. ■ 

We notice that the DoF region V coincides with the DoF region assuming perfect cooperation between the relays 
in each layer and there is no penalty due to distributed relays in terms of DoF. This property can be used to 
characterize the DoF region of more general networks having multi-antenna nodes and general message sets. Fig. 
[2] plots the DoF region V for K = 3 in which the sum DoF is limited by the minimum number of nodes in each 
layer. 
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Fig. 2. DoF region T> for K = 3. 
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Theorem 2: The DoF region V(W g ) of the Gaussian relay network is the set of all {dji}w oi ew g suc h that 



^ d^/72 < Lj :M +i for j E {1, • • • , i^M+i}, 

i=l 

Km+i 

^2 d J iI J i - Li ^ for i G {1, • • • , #1}, 
j=i 

V V djilji < min <^ V L^ m ) , 

^^ m€{l,-,M+l} f-' 



(13) 
(14) 
(15) 



where Iji is the indicator function that is one if Wji G W g and zero if W^ ^ W^. 

Proof: We refer to Section IV-BI for the proof. ■ 

The DoF region V(W g ) again coincides with DoF region assuming perfect cooperation between the relays in 

each layer. Fig. [3] plots V(W g ) when K x = K M+ i = 3, W g = {Wu, W 2 2, W33}, and L^ = L iiM +i = 2 for all 

i G {1, • • • ,3}. Note that this network is the same as that in Fig. [2] except that each node has multiple antennas. 

Specifically, the number of antennas at each source and destination is two. 



IV. ACHIEVABILITY FOR TWO-HOP AND THREE-HOP NETWORKS 

To prove the main results, we first study the two-hop and three-hop Gaussian relay network with K relays in 
each layer. We propose a block Markov AF relaying scheme and derive its achievable rate region, which will be 
used to characterize the DoF region in Section |V| 

A. Achievable Rate Region for M — 2 and K — K\ — K2 = K3 

In this subsection, we derive an achievable rate region of the Gaussian relay network for M = 2 and K — K\ — 
K 2 = K 3 . 

1) Channel space partitioning and pairing: As shown in the Introduction, interference-free communication is 
possible for every S-D pair if the relays amplify and forward their received signals with an appropriate delay r 
such that H.2[t + T ]Hi[t] becomes a diagonal matrix with non-zero diagonal elements. The relays, however, will 
have to wait forever in order to pair two channel matrices perfectly since channel coefficients vary according to 
a continuous distribution. To resolve this problem, we first partition channel spaces into subspaces and then pair 
subspaces of the first hop to those of the second hop. 

Define Q 4 {H 5 || real([H^-)| < 8Q, | imag([H^-)| < 8Q,iJ G {l,--- ,K},H 5 G 5(Z KxK + jZ KxK )}, 
where 6 > is the quantization interval and Q G Z + is related to the number of quantization points, i.e., card(Q) = 
(2Q + l) 2K . We further define %f as the set of all full-rank matrices in C KxK . Then we partition the channel 
space of the first hop into Hi(Hs) 9 s such that Hi(H 6 ) = {H| | real([H - U S ]ij)\ < S/2, | imag([H - H S ]ij)\ < 
5/2, i, j G {1, • • • , K}, H G Hf}, where H# G Q. To partition the channel space of the second hop and pair with 
%i(H$)'s, we define F : Hf — >• Hf such that 

F(II) ± c(H)II-\ (16) 

where c(H) > is set to satisfy p H2 [t](^( H )) = PH 1 [t]( H ), i.e., 

J] p^(c(H)|H- 1 ]y)= II Ph(\H\a)- ( 17 ) 

Then define « 2 (H$) = {F(H)|H G Hi(H s )}, which is the image of Hi(H s ) under F. 

Note that c(H) can be uniquely determined from H under the considered class of channel distribution. Also 
F(H)H becomes a diagonal matrix with non-zero diagonal elements and there exists one-to-one correspon- 
dence between H and F(H). Along with p Ha [ t ](F(H)) = p H i[t](H), this property will be used to prove that 
Pr(Hi[ti] G Hi(Hs)) = P r (H2[^2] G H2 (H^)) in Lemma Q] After the sources transmit their signals at time £ 
satisfying Hi [t] G ^(H^), the relays will amplify and forward them with the delay r such that H2 [t+r] G 7^2 (H#) 
(see Fig. @]). For this, we apply block Markov encoding and relaying. 

2) Block Markov AF relaying: Let us divide a length n block into B sub-blocks having length tlb = § each. 
We apply block Markov encoding and relaying over two hops so that the number of effective sub-blocks is equal 
to B — 1. Let 7m (H<$) be the set of time indices whose channel instances belong to % m (H$). That is, 

Ti b) (H d ) 4 {t|H m [t] G H m (H 5 ), t G {(6 - l)n B + 1, ■ ■ ■ , bn B }}, (18) 

where m G {1,2} and b G {1, • • • ,£>}. During the b^ effective sub-block, the sources and relays transmit using 
N(Hs) G Z + time indices in 7f (H$) and 7^ (H<$), respectively. The detailed procedure of the 6 th effective 
sub-block is as follows: 



• (Encoding) For all H# G Q, if card (7^ (H$)) < iV(H$) declare error, otherwise the sources transmit their 

T s) 



messages with power P using iV(H$) time indices in 7^ ^(H^ 



(Relaying) For all H5 G Q, if card (7^ (H«y)) < iV(H^) declare error, otherwise the relays amplify and 
forward their received signals that were received during 7f (H$) using iV(H$) time indices in 7^ (H5). 
Specifically, ^[h] = 7yi[^i] ? where t\ G 7f (H-s) and £2 G 7^ (H#). Here, the amplification factor 7 is 
given by 

7 ' 1+ ma* {II[HiM],II 2 f }p) ' <19) 

zG{l,---,A) / 




n 2 (u 5 ) 

Hi[t] H 2 [i + r] 



Fig. 4. Channel space partitioning and pairing. 



which satisfies the power constraint. 
• (Decoding) The destinations decode their messages based on iV(H$) received signals for all H^ G Q. 

Let E[ \, E^ \, and E^ \ denote the encoding, relaying, and decoding error events of the v^ S-D pair at the b^ 

effective sub-block, respectively. Notice that E[- occurs if card (7^ (H^)) < iV(H$) for any H^ G Q and E^- 

occurs if card(7^ (H^)) < N(H.§) for any H# G Q. From the union bound, 

P e , < J2 Pr (E[ b J ) + £ Pr (Eg ) + J2 Pr (Eg ) • (20) 

6=1 6=1 6=1 

3) Achievable rate region: We derive the achievable rate region of the proposed scheme. We first introduce the 
following lemmas, which will be used to prove the main theorem. 

Lemma 1: For all t l9 t 2 , and H s G Q, Pr (Hi[*i] G Hi(B. s )) is equal to Pr (H 2 [£ 2 ] G ^(H^)). 

Proof: Let us first show that F is bijective. Since F(H) is in Hf, it is enough to prove that F is injective, 
i.e., F(H a ) ^ F(H b ) if H a ^ H 6 for all H a G H f and H 6 G «/. Assume that F(H a ) = F(H 6 ), equivalently 
H« = ^}H b . Since F(H°) = F(H>) means n^{i,-^}^([ H 1^) = n^e{i,-,K}^([ H ^i) andp^a) ^ 
P/i(&) if \a\ 7^ \b\, we get H a = H 6 , which contradicts the condition that H a ^ H 6 . In conclusion, F is bijective. 
Then we obtain 



Pr (Hxfa] G Wi(H ff )) = / pH l[tl] (H)dH 



HGHi(H 5 ) 

PH a[ta] (F(H))dH 

^H 2 [t 2 ](H)dH 

HeH 2 (H 5 ) 

= Pr(H 2 [*2]GW2(H ff )), (21) 

where the second equality holds since £>Hi[ti](H) = Ph 2 [£ 2 ](^(H)) and the third equality holds since ^(H^) is 
the image of T-Li(Hs) under F and F is bijective. In conclusion, Lemma Q] holds. ■ 

Lemma 2 (Csiszdr and Korner): The probability that 



U B 



< e (22) 



for all H^ G Q is greater than 1 - card(Q)/(4n jB e 2 ), where m G {1, 2} and 6 G {1, • • • , B}. 

Proof: We refer to Lemma 2.12 in [20] for the proof. ■ 

Lemma 3: Suppose that H G ^(H^) and F(H) + A G ^(H^), where A is the quantization error matrix with 
respect to F(H). As the quantization interval 8 — >> 0, A converges to the zero matrix. 

Proof: From the definition of / H 2 (H 5 ), there exists H + £ G ^(H^) that satisfies F(H + S) = F(H) + A. 
Since |[53]^-| < 2\/25 — >> as 5 — » 0, SI converges to the zero matrix as 5 — >► 0. Hence, A also converges to the 
zero matrix as 8 — >► 0, which completes the proof. ■ 

Based on Lemmas Q] to [3l we derive the following theorem. 



Theorem 3: For the Gaussian relay network with M = 2 and K = K\ = K 2 = K3, 

^ = K los ( 1 + TTTOfFw )) <23) 

is achievable for all i G {1, • • • ,if}, where 7 2 = P/(l + max iG | 1 ... ^} {||[H]^||^} P) and c(H) is set to satisfy 
([T7I) . Here, the distribution of H is given by ]J i G r x ... K yPh([H]ij). 

Proof: Since Pr (H m [t] E % m (H$)) does not depend on t or m (Lemma [T]), we will use the shorthand notation 
Pr(H 5 ) to denote Pr(H m [t] E H m {R 5 )). From Lemmas CD and El we set N(H 6 ) = n B (Pr(H (5 ) - e). Then both 

Pr(E[ b J) and Pr(^) are upper bounded by ^ff , where we use card(Q) = (2Q + l) 2 ^ 2 . 

Suppose that the sources transmit their signals at t\ E 7f (H$) and the relays amplify and forward them at 
£2 £ T 2 (H<$). Let us denote Hi[ti] = H and H2[t2] = ^(H) + A, where A is the quantization error matrix 
with respect to F(H). From (|5]), 

y 2 [t 2 ] = 7 (c(H)I + AH) xi [h] + 7 (c(H)H- 1 + A)z x ft] + z 2 [t 2 ] , (24) 

where we use x 2 [£ 2 ] = 7Yi[ti] and F(H)H = c(H)I. Then the received signal-to-interference-and-noise ratio 
(SINR) at the z** 1 destination is given by 

SINR, = — ^KH) + [AH],,| 2 P (25) 

7 2 E l[AH]y|V+l + 7 2||[ c(H )H-i+A],|||. 

Hence, assuming that each source uses the Gaussian codebook, an achievable rate of the i^ S-D pair at the 6" 1 
effective sub-block is lower bounded by 

Rf ] > V min log(l + SINRi)(Pr(Ha)-e) (26) 



^ Q uen l{ u s ) 



with an arbitrarily small probability of decoding error, i.e., P[E^l) —¥ as tlb —> 00. Set 6 = n B '^ , 
Q = n^ , and e = n^ / , which are functions of n#. Then 

8 = n~ B ll{2AK2) -> (27) 

e = n~ 1/3 -> (28) 

5Q = nf 2 ^ -4 00 (29) 

(gQ + lg < 3^ n - 1 /a^ 

4n B e 2 " 4 " B ^° (3U) 

as n^ increases. The first condition guarantees an arbitrarily small quantization error, the second condition guarantees 
an arbitrarily small rate loss due to the randomness of channel realizations, the third condition is needed to use 
almost all channel instances for transmission, and the fourth condition guarantees arbitrarily small probabilities of 
encoding and decoding error. From (1251) to (l30l) and the result of Lemma [3] , we finally derive that 

lim iif > / log (l + - at/www-il 112 ) PH mW (H)dH (31) 

n B ^oo J HeHf \ 1 + 7 V(H)||[H ^iWjrJ N 

is achievable with probability approaching one. Since we can make both ub and B sufficiently large as n increases, 
Ri = ^jp^R\ —> R\ with P e ^ -^Oasn increases. From the fact that channel coefficients are i.i.d. drawn from 
a continuous distribution, the corresponding channel matrix is in %f with probability one. In conclusion, (1231) is 
achievable, which completes the proof. ■ 

The following example shows the achievable sum rate of the proposed scheme under i.i.d. Rayleigh fading. 

Example 1: Suppose the Gaussian relay network with M = 2, K = K\ = K2 = K3 and channel coefficients 
are i.i.d. drawn from A/*c(0, 1). Then, c(H) is given by ||H||i?/||H _1 ||i?. From the cut-set upper bound, one can 
derive a simple sum rate upper bound. For the cut dividing the sources and the rest of nodes, 

y^Ri<E\ max logdet(I + PHi[i]£ x [i]Ht[i])| , (32) 
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Fig. 5. Sum rate comparison for i.i.d. Rayleigh fading with M = 2 and K — K\ — K2 — K3. 



where £#[£] = ^pE(x[t]x^[t]) is the normalized covariance matrix. Fig. \5\ plots the achievable sum rate in Theorem 
[3] and its upper bound in (l32l) for K = 2, 3, 4. Numerical results suggest that, at any SNR, sum rate gaps are upper 
bounded by a constant independent of SNR. 



B. Achievable Rate Region for M = 3 and K = K \ = • • • = K4 

In this subsection, we derive an achievable rate region of the Gaussian relay network for M = 3 and K = 
K\ = • • • = K4. Since the basic concept and approach are similar to those in Section HV-A1 we mainly explain the 
proposed scheme and its achievable rate region and refer to the Appendix for the detailed proof. Together with the 
result of the two-hop network, the derived achievable rate region will be used to characterize the DoF region in 
Section [V] 

1 ) Channel space partitioning and pairing: Let T-Ld be the set of all full-rank matrices in C KxK whose 
eigenvalues are distinct. By using the eigenvalue decomposition, H e Hd can be represented as SAS -1 , where 
S consists of K eigenvectors and A = diag(Ai,-- , Ax) is the diagonal matrix with K ordered eigenvalues 
|Ai| > ••• > \\k\ > 0- Leti^ be the sign bit of real(A^), which is zero if real(A^) > and one if real (A^) < 0. Define 
(wi,w 2 , • • • , w K ) as (v K , v K vk-l ••• i v k @ vk-1 © * * * © ^i), A x = diag ((-l)™ 1 A x , • • • , (-1) Wk X k ), and 
A 2 = diag ((-l)™^ 2 , • • • , (-1)^A~ 2 ), where denotes XOR. Now define F x \U d ^ U d and F 2 :U d ^ %d 
as 

Fi(H) = ci^SAiS -1 and F 2 (H) = c 2 (H)SA 2 S- 1 J (33) 

where c m (H) > is set to satisfy Pu m+1 [t]( F m(H.)) = Ph^H), i.e., 

H P/l (c m (H)[SA m S- 1 ], J )= II PhWij)' 04) 

i,je{i,-,K} i,je{i,-,K} 

We partition the channel space of the first hop in the same manner as was done for HiCHs) and then partition and 
pair the channel spaces of the second and third hops based on F\ and F 2 . Specifically, we define 7^_(H$) = %i(H$), 
H' 2 (H S ) 4 {Fi(H)|H G H[(H S )}, and H' 3 {H 6 ) ± {F 2 (H)\H g H[(H s )}. 

Similar to F in Section HV-A[ c m (H) can be uniquely determined and F 2 (H)Fi(H)H becomes a diagonal matrix 
with non-zero diagonal elements. The role of (wi, • • • , wk) is to guarantee the one-to-one correspondence between 
A and A m , which also guarantees the one-to-one correspondence between H and F m (H). Without (wi, • • • , wk), 
2 K different matrices of A give a same diag (Aj -2 , • • • , A^ 2 ) (see Fig.[6]for K = 2). Along with Pu m+1 [t] (^m(H)) = 
PHi[t](H), this property will be used to prove that the probabilities of paired channel subspaces are the same in 
Lemma HI 
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Ai A 2 

-• 



diag(a, b) diag(a, b) diag(a 2 ,6 



2 h-2\ 



-• • 



diag(a, -b) diag(-a, b) diag(-a 2 ,-6 2 

• • • 

diag(— a, b) diag(— a, — b) diag(a -2 , — b~ 2 ) 

• • • 

diag(-a, -b) diag(a, -b) diag(-a~ 2 , b~ 2 ) 

A diag(Ai,A 2 ) diag(A^ 2 , A^ 2 ) 

• • ^# 

diag(a, b) diag(a, b) ^^diag(a~ 2 , b~ 2 ) 



diag(a, -b) diag(a, -by /diag(a 2 ,-6 2 ) 

• v / • 

diag(-a, b) diag(-a, b) diag(-a" 2 , b~ 2 ) 

• m • 

diag(— a, — b) diag(— a, — b) diag(-a~ 2 , — b~ 2 ) 

Fig. 6. The mapping functions Ai and A2 for K — 2 (upper) and the mapping functions diag(Ai, A2) and diag(A^~ 2 , A^ 2 ) for K — 2 
(lower), where real (a) > and real (b) > 0. 

2) Block Markov AF relaying: We divide a length n block into B sub-blocks having length ub = § and then 
apply block Markov encoding and relaying. Thus, the number of effective sub-blocks is B — 2. Define 

7^(Ha) 4 {t|H m [t] E ?C(Htf),t G {(6 - l)n B + 1, ■ ■ ■ ,6n B }} , (35) 

where m E {1, 2, 3} and 6 E {1, • • • , B}. During the b^ effective sub-block, the sources, the relays in the second 
layer, and the relays in the third layer transmit using N'(Hs) E Z + time indices in T[ (H$), 7^ (H«j), and 
7^ (H<y), respectively. The detailed procedure of the 6 th effective sub-block is as follows: 

• (Encoding) For all H5 E Q, if card (7^ (H«$)) < N'(Hs) declare error, otherwise the sources transmit their 
messages with power P using N'(Hs) time indices in T[ (H$). 

• (Relaying) For all m E {1,2} and H 6 E Q, if card(7^ m) (H<5)) < N f (H 6 ) declare error, otherwise the 
relays amplify and forward their received signals that were received during 7™ m ~ (H$) using N'(Hs) time 
indices in 7^ m) ( H <0- Specifically, x m+ i[£ m+ i] = J m y m [t m ], where t r E 7^ /(6) (H 5 ), t 2 E 7^ (6+1) (H,5), and 
£3 E 7^ (H5). Here, the power amplification factors are given by 

/ P v /2 

71 ll+ max A\\mhU 2 F }p) ' } 

and 

72= I 1 + 7?. max ^{||[H 2 [t 2 ]] i |||, + ||[H 2 [t 2 ]H 1 [t 1 ]] i |||.P} J ' (3?) 

\ zG{l, •••,#} / 

which satisfy the power constraint. 

• (Decoding) The destinations decode their messages based on N f (Hs) received signals for all H# E Q. 

3) Achievable rate region: We derive the achievable rate region of the proposed scheme. Because the overall 
proof is similar to that of Theorem [3l we explain the main result here and refer to the Appendix for the detailed 
proof. 
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Fig. 7. Sum rate comparison for i.i.d. Rayleigh fading with M — 3 and K — K\ 



K 4 . 



Theorem 4: For the Gaussian relay network with M = 3 and K = K\ = • • • = K4, 

7 2 2 7i 2 ci(H)cKH)P 



i*i = £ log 1 + 



l + Tlll[^2(H)],||2, + T 2 T 2|| [jP2(H)j p l(H)] 



illF 



is achievable for all i G {1, 

w2 



,/f}, where 7^ and 7I are given by P/(l + max^^ 

12 1 lirz? /-tJMJl.112 ^ 



P/(l + i( max, 6{1) ...^ } {||[Fi(H)]i|||, + ||[Fi(H)Hy2,P}), respectively. The definition of F m (H) is given by 



rf . 

,K}P/i([ H Ju 



H ]dlF 



(38) 



}P) and 



(l33l) and c m (H) is set to satisfy (l34l) . Here, the distribution of H is given by Yl i G r x 

Proof: We refer to the Appendix for the proof. ■ 

Example 2: Suppose the Gaussian relay network with M = 3 and K = K\ = • • • = K4. If channel coefficients 
are i.i.d. drawn from A/c(0, 1), then ci(H) and C2(H) are given by |H||i^/|| SAiS - - 1 ||,p and |H||i?/||SA2S _1 ||i?, 
respectively. Fig. [7] plots the achievable sum rate in Theorem |4] and its upper bound in (l32l) for K = 2, 3, 4. Although 
the sum rate gaps are larger than those of the two-hop network, numerical results still suggest that, at any SNR, 
sum rate gaps are upper bounded by a constant independent of SNR. 



V. DoF Region 

In this section, we characterize the DoF region V of the Gaussian relay network based on the results in Section 
HVl We further characterize the DoF region V(W g ) assuming multi-antenna nodes and general message sets. 



A. DoF Region of Gaussian Relay Networks 

In this subsection, we prove Theorem [TJ The converse can be shown from a simple cut set upper bound. From 
the cut dividing the z tJl source and the rest of nodes, Ri is upper bounded by 1 x K 2 single-input single-output 
(SIMO) capacity, which gives d{ < 1, where i G {!,••■ , K}. From the cut dividing the nodes in the m tn layer 



\th 



^K 



and the nodes in the (m + 1) U1 layer, ^ =1 R% is upper bounded by K m x K m +i multiple-input multiple-output 
(MIMO) capacity, which gives J2i=i &i ^ min {^m 5 K m +i}, where m G {1, • • • , M}. Hence we obtain (fTTT) and 

Let us now show the achievability. We first prove that every S-D pair achieves one DoF if K = K\ = • • • = 
Km+i- Consider the case where M = 2. We recall the rate expression (l23l) in Theorem [3] that is lower bounded by 

7 2 c 2 (H)P 



R z > E log 1 + 



l + 7 2 c 2 (H)||H- 



(39) 



Ri > E ( log ( 1 + , , : ,„ - ^ 1 '\,'l' 2 I, E ^,,2 ) ) • (4<» 
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Assuming P > 1, we derive 1/(1 + ||H||^) < 7 2 < if/||H||^, showing that the bounds are independent of 
P. Because the absolute values of channel coefficients are non-zero and finite with probability one, the re- 
lated quantities such as c 2 (H), ||H||^, and ||H _1 ||^ are also non-zero and finite with probability one. Hence 
limsupp^QQ-Ri/logP = 1. Since the above bounds hold for any z, every S-D pair achieves one DoF. Consider 
the case where M = 3. We recall the rate expression (l38l) in Theorem |4] that is lower bounded by 

72 2 7l 2 cj(H)c 2 (H)P 

i + tIII^^iii + t^ii^^iiiiip^h)!!^ 

Assuming P > 1, we derive 1/(1 + ||H||2,) < 7 2 < K/\\U\\ 2 F and 1/(1 + 7i(||Fi(H)|||. + ||Fi(H)|||.||H|||)) < 
72 < ^/(7i ll^i(H)H||p). Note that the bounds are again independent of P. Hence limsupp^^ Ri/logP = 1. 
For M > 4, one can immediately obtain that every S-D pair achieves one DoF from the results of M = 2 and 
M = 3. If M = 5, for example, we can apply the result of M = 2 over the first two hops and the result of M = 3 
over the rest three hops. 

Let us now consider the case where K m ^ K in general. In this case, V has corner points (d*, • • • , d* K ) such 
that Ya=i d k = min m {K m } and d\ e {0, 1} for all i G {1, • • • , if}. To achieve (d|, • • • , d* K ), min m {if m } S-D 
pairs with d* = 1 participate in communication and min m {if m } relays in each layer participate in relaying. Hence 
one DoF is achievable for each of the corresponding min m {if m } S-D pairs, where we use the result for the case 
K = K i = • • • = Km+i- Note that any point in the dominant face can be achieved by time sharing between corner 
points. In conclusion, Theorem Q] holds. 

B. Multi-antenna and General Message Set 

In this subsection, we prove Theorem [2 The converse can be shown from the cut-set upper bound. From the 
cut dividing the j^ destination and the rest of nodes, J2i=i Rjilji * s u PP er bounded by (X^Ji ^i,m) x £j,m+i 
MIMO capacity, which gives (fT3l) . From the cut dividing the v source and the rest of nodes, J2j=i +1 Rjilji * s u PP er 
bounded by L^\ x {J2j=i Ljp) MIMO capacity, which gives (fT4l) . Lastly, from the cut dividing the nodes in the rrv 

layer and the nodes in the (ra+l) th layer, Ya=i J2j=i +1 R jilji is upper bounded by (J^i L^ m ) x (^ • J^ +1 L^ m+i ) 
MIMO capacity, which gives (IT5l) . 

Consider the achievability. The greedy allocation of {dji}w Jt ew g based on a specific order that satisfies (fT3l) to 
([T5l) is one of the corner points in V(W g ). Let {rf*J^.. e yy g be the result of the greedy allocation according to a 

specific order. For Wji G )%, we can choose d*^ antennas at the i^ source and d*^ antennas at the j th destination 
and pair them as d*^ virtual S-D pairs. As a result, we can establish total X^=i S 7 ^i +1 ^ji^ji virtual S-D pairs 
because {d^ i }w 3% ^w g satisfies (IT3T) and (fT4l) . We can also choose total J2i=i J2j=i +1 d^Iji relay antennas in each 
layer because {d^jw^eWg satisfies ([T5l) . The resulting network consists of J2i=i J2j=i +1 d%Iji virtual S-D pairs 
with J2i=i J2j=i +1 djjlji virtual relays in each layer. Therefore, from the result of Theorem [TJ all virtual S-D pairs 
can achieve one DoF, equivalently {d^ j }w 3% ^w g * a achievable. Note that any point in the dominant face can be 
achieved by time sharing between corner points, which completes the proof. 

VI. Conclusion 

In this paper, we study the if -user M-hop Gaussian relay network consisting of K m nodes in the rmr 1 layer. 
The proposed block Markov AF relaying exploits channel fluctuation to cancel the inter-user interference at each 
destination and works for a wide class of channel distributions including i.i.d. Rayleigh fading. Under this class 
of channel distributions, we completely characterize the DoF region of the Gaussian relay network and further 
characterize the DoF region of more general networks with multi-antenna nodes and general message sets. The 
resulting DoF regions coincide with the DoF regions assuming perfect cooperation between relays in each layer 
and there is no penalty due to distributed relaying in terms of DoF. 

Appendix 
Achievable Rate Region for M = 3 and K = K\ = • • • = K 4 

In this appendix, we derive the achievable rate region in Theorem 01 The overall steps are similar to those 
in Section IIV-AI We first introduce a useful property, which will be used to prove the following lemmas. Let 
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H a = s a A a (S a )- 1 and H 6 = S 6 A 6 (S 6 )-\ where H a G U d and H 6 G H d . Then H a = H 6 if and only if 
A a = A 6 and there exists a diagonal matrix satisfying S a = S 6 0. Similar to Lemmas [T] and [3 we derive the 
following two lemmas. 

Lemma 4: For all ti, £2, ^3? an< 3 H# G Q, 

Pr(Hi[ti] G «i(H*)) = Pr(H 2 [t 2 ] G «£(H*)) = Pr(H 3 [* 3 ] £ n' 3 (H s )). (41) 

Proof: The overall proof is similar to that of Lemma [TJ To show that F m is bijective, we need to prove 

F m (H a ) ^ F m (H 6 ) if H a ^ H 6 for all H a G « d and H 6 G « d . Assume that F m (H a ) = F m (H b ). Then 

S a A^(S a ) _1 = 4^JS 6 A^(S 6 ) -1 , where S a and A a (S 6 and A 6 ) are the eigenvector matrix and the eigenvalue 

matrix of H a (H 6 ), respectively. Hence F m (H a ) = F m (H 6 ) if and only if A^ = 4j4a^ and there exists a 
diagonal matrix satisfying S a = S 6 0. Since there is one-to-one correspondence between A and A m , we get 
A a = cA b , where c > is a constant. As a result, we obtain 

H a = S a A a (S a )- 1 = cS^A^S 6 ©)- 1 = cH 6 . (42) 

Since F m (H a ) = F m (H 6 ) means Uije{i,- ,K}Ph(\H a ]ij) = U i9J e{i^ 9 K}Ph(\^ b ]ii) ^d p h {a) ^ p h (b) if \a\ + 
|b|, we get H a = H 6 , which contradicts the condition that H a 7^ H 6 . In conclusion, F is bijective. Then similar 
to (|2Ti we obtain Pr(Hi[*i] G ^(H^)) = Pr(H m [£ m ] G ^(H^)), which completes the proof. ■ 

Lemma 5: Suppose that H G 7^i(H 5 ) and F(H) + A m _i G ^(H^), where m G {2,3}. As 5 -» 0, A x and 
A 2 converge to the zero matrix. 

Proof: Since the overall proof is the same as that in Lemma we refer to the proof of Lemma [51 ■ 

We are now ready to prove Theorem 0] by using the previous lemmas. The overall procedure is similar to that of 
Theorem [3j Since Pr (H m [t] G T-L' m {H.s)) does not depend on t or m (Lemma |4]), we will again use the shorthand 
notation Pr(H 5 ) to denote Pr(H m [i] G H' m (B.s)). We set N'(H) = n B (Pr(H 6 ) - e). Then the probability of 

th 3(2Q-\-1) 2k2 

encoding and relaying error of the b in effective sub-block is upper bounded by v ^~ e ; 2 — . 

Suppose that the sources transmit their signals at t\ G T[ (H$) and the relays amplify and forward them 
at t 2 G 7^ (6+1) (H 5 ) for the second hop and £3 G T^ b+2 \n 6 ) for the third hop. Let us denote Hi[*i] = H, 
H2[t2] = Fi(H) + Ai, H^^] = F 2 (H) + A2, where Ai and A2 are the quantization error matrices with respect 
to -Fi(H) and F 2 (H), respectively. Then 

ys[ts] = 7271 (c 2 (H)ci(H)I + A tot ) xi[*i] 

+ 7271 (F 2 (H) + A 2 )(Fi(H) + Ai)zi[ti] + 72 {F 2 (H) + A 2 )z 2 [t 2 ] + z 3 [t 3 ], (43) 

where 

A« = F 2 (H)AiH + A 2 F!(H)H + AsAxH. (44) 

Hence the received SINR at the i^ destination is given by 

SINK, _ ^»f Ma) + '^l' F (45) 

where a\ oU = 1 + 7 2 2 H[^(H) + A 2 ]i\\ 2 F + 7 2 2 7 2 ||[(^(H) + A 2 )(F 1 (H) + Ai)]^*.. Hence, 

i?f ) > V min log(l + SINR i )(Pr(H 5 )-e) (46) 

^ H6«J(H,) 
is achievable with an arbitrary small probability of decoding error. Setting 5 = n^ ^ \ Q = n^ , and 



n 



-1/3 



£ 



satisfies the conditions in (l27l) to (l30l) . Hence, 

lim i? W 

TIB— ^OO 

is achievable with probability approaching one, where we use the result of Lemma [5J Since we can make both tlb 
and B sufficiently large as n increases, Ri = ^j^R^ ^ R\ with P ei -^Oasn increases. In conclusion, 
is achievable, which completes the proof. 
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